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ANALYSIS OF EFFECTS OF SPANWISE VARIATIONS OF GUST VELOCITY 
ON A VANE-CONTROLLED GUST- ALLEVIATION SYSTEM 

By L. Keith Barker and Gene W. Sparrow 
Langley Research Center 

SUMMARY 

An analysis has been made of the effects of spanwise variations of gust velocity in 
isotropic turbulence on a gust- alleviation system which employs an angle-of-attack vane 
mounted ahead of the wing to sense the vertical gust velocity. The wing flaps were moved 
in response to the vane deflection by a linear second-order servosystem to produce a lift 
opposite to that produced by the gust. Since the gust velocity actually varies across the 
wing span and the vane only senses the gust velocity at a single point, the flaps cannot be 
expected to eliminate the effects of the gusts completely by being deflected in response to 
a signal from the vane. The effectiveness of the alleviation system is examined by con- 
sidering the statistical properties of the residual lift on the alleviated wing. Considera- 
tion of spanwise variation of gusts has indicated design parameters (gain and natural fre- 
quency of the flap servosystem and vane location) that give substantial reduction of the 
lift due to gusts. 


INTRODUCTION 

To insure passenger comfort and reduce structural loads when an airplane flies 
through rough air, it is desirable to minimize resulting accelerations. In 1951, a rather 
thorough theoretical analysis of some methods for increasing the smoothness of flight 
through rough air was made. (See ref. 1.) As a result of this analysis, a method of gust 
alleviation by direct lift control evolved and was later flight tested. (See refs. 2 and 3.) 
Although the gust alleviation was substantial in the flight tests, the system did not work as 
well as had been anticipated. 

With the increasing emphasis in space projects, this gust alleviation project was 
discontinued. Also, other reasons why there was less interest in gust alleviation at this 
time are (1) jet transports tended to fly at high altitudes which reduced the frequency of 
gust encounters, (2) higher wing loading was used, and (3) improved weather forecasting 
and the ability to locate storm centers permitted aircraft to avoid turbulent air. 



Interest in gust-load alleviation has recently begun to increase because of possible 
applications to V/STOL aircraft. This in turn has called attention to the flight tests men- 
tioned previously. 

There are several possible reasons why less gust- load alleviation was obtained in 
the flight tests (refs. 2 and 3) than had been predicted in reference 1. Three of these are 
(1) the effect of spanwise gust variation, (2) inadequate or nonlinear servo response, and 
(3) nonlinear flap effectiveness. In the present study the reduction in the wing lift due to 
turbulence is examined with the spanwise variation of the gusts taken into account. Some 
results are also included of the desired characteristics of the flap servomechanism. 

SYMBOLS 


Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. 

b wing span 


wing chord 

wing mean aerodynamic chord 

nondimensional gust frequency, 

lift coefficient 

local lift coefficient 

8Cl 


C 

Cl 


-‘OL 


‘'a 

Cl 


6f 


CjjgL 


lift- curve slope. 
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gain constant, — — =■ Kj 

llQ, 

gain constant in flap control system 

COr|L 

nondimensional natural frequency of flap servo, ” 
scale of turbulence 

vane location forward of aerodynamic center of the wing 

nondimensional vane location, 

lift on wing due to flap deflection 

lift on wing due directly to vertical gusts 

dynamic pressure 

percentage reduction in the power spectrum of the basic wing when the 

<p-t ~ 

alleviation system is used, — ~ — — x 100 

01 

percentage reduction in mean- square lift on basic wing when the alleviation 

01 ( 0 ) - 0 .( 0 ) 

system is used, ^ — x 100 

wing surface area 
time 

constant forward velocity of the wing 
vertical gust velocity 
reference axes 
coordinates 

y 

nondimensional spanwise coordinate, 

b/2 
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angle of attack on the wing 


angle of attack due to vertical gust velocity 
angle of attack of vane due to vertical gust 
lift weighting function 
Dirac delta function 
flap deflection 

damping ratio in flap control system 
dummy variable of integration 

Ut 

nondimensional distance along flight direction, — 

Xj 

distance between two points in space 

time increment for wing traveling at constant velocity 
from Xj to X2 

power spectrum of e(x) 

nondimensional power spectrum of e(x) 

nondimensional power spectrum of 

correlation function of e(x) 

nondimensional correlation function of e(x) 

correlation function of vertical gust 

nondimensional correlation function of 


to travel 


spatial natural frequency of flap servo in radians per unit of distance 





gust frequency 


ci)n natural frequency of flap servo 

All primed values denote nondimensionalization by the scale of turbulence L. For 
example, b' = ^, The symbol A preceding a quantity means an incremental value. 

Xj 

For example, Ax = X 2 - xj. 


STATEMENT OF PROBLEM 

The problem being examined is illustrated in figure 1 . The airplane is assumed to 
be flying in a straight line with wing level through gusty air which varies not only longi- 
tudinally but also across the wing span of the airplane. Based upon the gust sensed at a 
point by a small vane, the flaps are deflected either up or down as required to alleviate 
the effect of the gust. In some cases, the flaps were assumed to deflect in phase with the 
vane, and in others, the flaps were assumed to be operated by a second-order servomech- 
anism. Note that the gust velocity hitting the vane is not necessarily equal to the gust 
velocities at various points along the span. For this reason, the flap system cannot be 
expected to eliminate the gust load totally if it responds to the vane deflection. A mea- 
sure of the residual lift on the wing can, however, be estimated by using the statistical 
properties of gusts to calculate the statistical properties of the residual lift. 

A method for computing the statistical properties of lift on an airplane wing due to 
random atmospheric turbulence is given in reference 4. The crux of the approach is to 
consider the turbulence isotropic so that the statistical characteristics of the turbulence 
are independent of direction. In this manner, the statistical properties of the vertical 



Figure 1.- Airplane flying through isotropic turbulence. 
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gust velocities along the wing span can be related to those of the vertical gust velocities 
along the flight path, and consequently, an analytical solution to the problem can be 
formulated. 

The approach of reference 4 is used in the present study to examine the statistical 
properties of lift on a wing due to vertical gusts when the wing is equipped with an allevia- 
tion system. 


ANALYSIS 

The coordinate system used in this study is shown in figure 2. A rigid wing is 
moving along the X-direction with constant forward velocity U, A vane is mounted a 
distance I ahead of the aerodynamic center to sense the gust. 



Direction of flight 

t 
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^Wing 
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pVane sensor 




> 

1 
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I 

1 

1 

1 
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A ^ 1 


Y 

Figure 2.- Coordinate system used In this study. 


Lift on the Wing- Flap System 

A theorem in linearized aerodynamics states that the lift in steady or indicial 
motion of a wing having arbitrary twist and camber is equal to the integral over the plan- 
form of the product of the local angle of attack and the loading per unit angle of attack at 
the corresponding point of a flat-plate wing of identical planform in reverse flow (ref. 5, 
pp. 325-326). The wing lift, therefore, can be written as 




unit a 


Q!(x,y*)dy* 


( 1 ) 


The local lift coefficient per unit angle of attack of the wing is equal to the local section 
lift-curve slope; therefore, equation (1) can be written as 
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( 2 ) 


L^(x) = ^ I j* ^ a(x,y*)dy* 

Multiplying and dividing through by the product of wing lift- curve slope and wing area 
results in 


qSCx pi ccj 

Lw(x) = — Z0-^Q!(x,y*)ciy* (3) 

Since the local section lift- curve slope and the total wing lift- curve slope are based on 
the same angle of attack, equation (3) can be written as 



1 cc, 

-^Q!(x,y*)dy* 
■1 cCl 

(4) 


1 

y(y*)a(x,y’^dy* 

.1 

(5) 


The angle of attack at location (x,y) caused by a vertical gust is given approximately by 

w(x,y) 


aJx,y) = 


U 


( 6 ) 


The wing lift caused by a vertical gust which varies across the span is therefore given by 

(7) 


qSCL . pi 

L^(x) = ^ ^ y(y*)w(x,y*)dy’ 


If the wing is provided with a flap, the increment in lift caused by flap deflection 
can be expressed by 


Lf(x) = qSCLg^6f(x) (8) 

In this analysis it is assumed that the flap deflection 6f(x) is the response of a 
linear system, with angle of attack of the sensing vane q!(x + Z, 0) as the input variable; 
therefore, the flap deflection is given by the familiar Duhamel convolution integral as 


6f(x) = Ki \ h(X)o! v(x + 1-X, 0)dX (9) 

-'0 

where h(x) is the flap response due to a unit impulsive vane angle of attack, and h(X) = 0 
for X < 0. In the section "Application of Equations," expressions will be introduced for 
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h(x) which correspond to the following cases: (1) the flap moving in phase with the vane 
deflection and (2) the flap deflection related to the vane deflection through a linear system 
of second order. 


By assuming no lag in the response of the vane to a gust, the angle of attack of the 
vane can be expressed as 

w(x + l, 0) 


Q!v(x + Z,0) = 


U 


( 10 ) 


Substituting equations (9) and (10) into equation (8) results in 

U Jo 


Lf(x) = Kj — \ h(x)w(x + i- A, 0)dX 


( 11 ) 


In the gust- alleviation concept of this paper, the flap is used to reduce the lift 
caused by the gust. The lift on the alleviated wing due to vertical gusts therefore is 
given by 

e(x) = L^(x) - Lf(x) 
or, from equations (7) and (11), 


e(x) = 

U 


I- r y(yV(x,y*)dy* - K f h(x)w(x + z- x,0)dx 
2J_i Jo 


( 12 ) 


Statistical Properties of Lift on a Wing Due to Vertical Gusts 

The statistical properties of interest are the correlation function and power spec- 
trum of the lift on the wing-flap system . These functions depend somewhat on the assumed 
model of atmospheric turbulence. In this study, it is assumed that the turbulence is sta- 
tionary, homogeneous, and isotropic. Under these assumptions the correlation function 
between vertical gust velocities at two points depends only on the distance between the 
points, and not on their location in space. The correlation function, therefore, can be 
expressed by 

1 

’^w{^2-^l»y2-yi) ^ W(xi,yi)w(x 2 ,y 2 )dxi 


= W^g 


\/(Ax) 2 + (Ay) 2 


= w2g(or) 


(13) 
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where w2 is the mean-square vertical gust velocity, and g(a) is a nondimensional 
correlation function for vertical gusts which is a scalar function of the distance between 
the two points (ref. 6). 

Correlation function of lift .- The correlation function for the resultant lift on the 
wing at points xj and %2 along the direction of flight is defined by 


where X 2 = Xj + Ax. Substitution of equation (12) into equation (14) and then making use 
of the relation shown in equation (13) results in 


'#'ep2-=‘i) =(%^) ”'^4 j"*! I.*i >-(yi*My2*)e +(|f 


hr ■1^2* 




dXg dyj 




dXi dyg’ 


k 2 h(Xi)h(X2)g(|Ax- Ax|)dXi dX2| 


(15) 


It is assumed that correlations in space and time are equivalent so that X 2 - xj 
can be replaced by Ut in equation (15). This assumption is a result of Taylor's hypoth- 
esis. (See ref. 6.) For later convenience this substitution is made and equation (15) is 
written in the following form: 

%{Vr) = - Kxly^iVr) - Kxjy^iVT) + (16) 

where 
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(18) 


1/^3 (Ut) 






ayj dyj 


'^rilo >'(>'1 * \j('^^*i- >^2? ■►(if (y 1 •/ 


dX2 dyj’ 




dXi dyg 


= j|) Jq h(Ai)h(X2)g(|UT-Ax|)dXidX2 


(19) 

( 20 ) 


( 21 ) 


Mean-square value of lift .- The mean-square value of lift is obtained by averaging 
the square of the lift e(x) throughout the region of turbulence, which is precisely the 
value obtained by evaluating equation (14) with Xj = X2> that is, ipjfi)- The nondimen- 
sional mean-square value of e(x), therefore, can be written as (from eq. (16)) 

5^e(0) = >^1(0) - k2,//^(0) (22) 

From the definitions of ^'^d \}/^ given by equations (19) and (20), it follows 

that when Ax (or Ut) is equal to zero. 


^2^0) = i//g(0) 

and equation (22) reduces to 


5^g(0) = )//j(0) - 2KiJy^i0) + 


(23) 


The term represents the nondimensional mean-square lift on the wing with 

no alleviation. The remaining terms are associated with alleviation obtained by flap 
deflection. Since the purpose of a gust- alleviation system is to make the total wing lift 
as small as possible, it is desirable to select the constant K to minimize i^g(O). The 
proper value of K can be obtained by differentiating equation (23) with respect to K, 
setting the result equal to zero, and solving for K. The result is 


K = 


1/^2 (®) 


( 24 ) 


Power spectrum of lift.- The power spectrum of wing lift e(x) is defined as the 
Fourier transform of the correlation function of e(x); that is. 
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Since the function xI/q{Vt) is an even function, equation (25) can be written as 

M^S) = ?io ^26) 


or, from equation (17), 


^Cx qS\" 


B = Kir"/ ^lo 


(27) 


where 


Wo-L 

<= = -5- 




Inspection of equation (16) shows that equation (27) can be written as the sum of 
four terms: 


where 


^e(C) = 0i(C) - Kcf>2iC) - K<^3(C) + K^(f>^{C) 


0e(C) = 


*e(c g) 


and 


e'-' - , ,2 

A \J u 


0i(C) = J°° i^.(L^)cos(C^)d^ 


(28) 

(29) 


(i = 1, 2, 3, or 4) 


(30) 


APPLICATION OF EQUATIONS 

In order to compute the correlation function, mean-square values, or power spec- 
trum of lift for a specific configuration, expressions must be assumed for the lift weighting 
function y(y), the function g(a) associated with isotropic turbulence, and the impulsive 
flap response h(x). 

The expression for y(y*) used in this study corresponds to an elliptically loaded 
wing and is given by (ref. 7) 
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(31) 


y(y*) = -/i - (y*)2 

A form of the function g(a), which has previously been considered in reference 7 
in studying isotropic turbulence, is used in this study. For two vertical gust velocities 
located at the points and ^X 2 ,y 2 ^, this function takes the form 

g(CT) = exp^-2 g-j (32) 

where 

^ = /(X2 - xi)2 +(y2 - Yif = ^(Ut) 2 +(|)^(y2* " yi*f (33) 

The flap response h(X) to an impulse depends on the assumed characteristics of 
the flap and its drive system and on its relation to the angle sensed by the gust- sensing 
vane. Equations will be derived for two types of flap response: (1) Instantaneous 
(in-phase) response and (2) response described by a second-order linear differential 
equation. 


In-Phase Flap Response 

When the flaps move in phase with the vane, h(X) becomes a Dirac delta function 
which is written as 


h(X) = 6(X) (34) 

The correlation function of e(x) is determined by use of equations (31) to (34) and equa- 
tions (16) to (21). The results are 

- Ki//3(^) + (35) 


where 




’l 

U / 



(36) 
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i// 4(0 = exp(-2|^|) (40) 

Nondimensional mean- square value of lift .- The nondlmensional mean-square value 
of lift is given by equation (23), which reduces to the form 

^^(0) = rpiiO) - 2Kily^{0) + (41) 

where 

- (y2*f - yi*)|'‘yr 

and 

Note that the nondimensional mean-square value is determined by the parameters b' 
and l\ 

With the substitution of i// 2 ( 0 ) and 1 ^ 4 ( 0 ) from equations (38) and (40), respec- 
tively, into equation (24) for the gain constant K which minimizes the mean-square value 
of e(x), the following result is obtained: 



Nondimensional power spectrum of lift .- The nondimensional power spectrum of 
lift is computed by using equations (27), (29), and (30). 

Effect of Lag in Flap Response 

In this study, it is assumed that the flap deflection 6^{x) is related to the angle of 
attack of the gust sensor o;y(x + f, 0 ) by means of a linear second-order differential equa- 
tion of the form 

2 

6f (x) + 2?n„ 6f (x) + (x) = KiS2„2q;^(x + z, 0) (45) 

where ^ is the damping ratio, is the natural frequency in radians per distance, and 
Kj is a gain or proportionality constant. 
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For a damping ratio ? < 1, the solution of equation (45) is given by equation (9) in 
which h(X) is equal to 


h(X) = -=S=exp(-COn^)sin(nn/l - ?2 xj 

/l - 


Since x = Ut, can be replaced in equation (46) by the equivalent expression 


fin = Wi- 


With this substitution, equation (46) becomes 




Correlation fimction of lift.- The correlation function of lift can be determined by 
substituting equations (31), (32), (33), and (47) into equations (16) to (21), and the result is 

= i'lU) - Ki// 2(4) - Ki//g(4) + k2^^(4) (48) 


= 


. u / 
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i) \/^ ■ (^2*f exp(-?kXi')sin(k/l - ?2 Xj’) 


dXj^' dyg* 


' tfri) ^0 io * >-2')]M'^f^ h) 

sin(k^l - ?2 X 2 ')exp^- 2 j^ + Xj' - X 2 'j^dX]^' dX 2 * 


(53) 


Mean-square value of lift .- The mean- square value of lift can be determined by 
evaluating equations (48) to (53) with ^ equal to zero. 

Power spectrum of lift .- The power spectrum of lift is determined by substituting 
equations (50) to (53) into equations (28), (29), and (30) and integrating the resulting 
expressions. 


RESULTS AND DISCUSSION 

A basic criterion used in this study to evaluate the effectiveness of the alleviation 
system is the reduction in the mean-square value of wing lift e(x). The mean square of 
e(x) is discussed, then the power spectrum of e(x) is used to examine how this reduc- 
tion in mean-square lift is distributed over the different gust frequencies. 

The various integrals in this study were evaluated numerically on a digital 
computer . 


In-Phase Flap Response 

Nondimensional mean- square lift .- Figure 3 shows the nondimensional mean-square 
lift AO) due to gusts acting on a wing flying through isotropic turbulence with and with- 
out the gust alleviation system. The nondimensional mean-square lift is presented as a 
function of b', the ratio of wing span to the scale of turbulence, and I*, the vane loca- 
tion in fractions of wing span. 

Note that the nondimensional mean-square lift on the basic wing (no alleviation) 
is reduced as the wing span becomes larger relative to the scale of turbulence 
(b’ increases). This self-alleviation is due to a spanwise averaging effect; that is, 
the lift increments due to various gust velocities tend to average out. 
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Figrire J.- Effect of alleviation system on the mean-sq.uare value 
of e(x). (Flap moving in phase with vane and optimum value 
of gain constant K ixsed.) 

The solid lines represent the results for the flaps moving in phase with the vane, 
that is, when there is no lag in the system. The vane senses the gust and the flaps 
respond immediately. The values of K used to obtain these curves were those which 
minimized the mean-square lift, as determined from equation (44). 

A wide range of b' values are shown in the figure to clarify the various trends; 
however, practical values of b’ appear to be rather small - on the order of 0.1 or less. 
For the flight tests in references 2 and 3, the vane location I* appears to be about 0.5. 
At this value of I*, 4/^(0) is reduced from 1 to 0, or 100 percent, at b' = 0 and from 
0.94 to 0.14, or 85 percent, at b' = 0.1. Thus, for practical values of b' up to 0.1, the 
figure shows that the nondimensional mean- square lift should be reduced by 85 to 100 per- 
cent, provided the response of the servo driving the flaps is fast enough. 

Gain constant K .- The results of figure 3 were generated by choosing a value of 
K to minimize the mean-square value of lift on the alleviated wing. The variation of gain 
constant with b' and I* is shown in figure 4 . Note that as the wing span and vane 
location ahead of the wing increase, K decreases. This would be expected since the 
gust velocity measured by a vane far ahead of the wing would not be expected to be useful 
in predicting the present gust velocity on the wing; hence, the optimum value of K would 
approach zero as I* approaches infinity. Likewise, as the wing span increases, the 
gust velocities on the wing could differ appreciably from the gust velocity measured by 
the vane at a point, and the optimum value of K would again approach zero; however, it 
turns out that the wing is alleviating itself in this case as shown previously in figure 3. 

For practical values of b’, the gain constant varies from about 0.8 to 1.0. 


16 




Figure 4.- Gain constant used to minimize the mean-square value 
of e(x). (Flap moving in phase with vane.) 

Power spectrum of lift .- It is of interest to examine the reduction in the power 
spectrum over the gust frequency range which results from minimizing the mean square 
of e(x). The nondimensional power spectrum of the lift on the basic wing is plotted 
against the nondimensional gust frequency C on a log- log scale in figure 5 for b' val- 
ues of 0.05, 0.10, and 1.0. As expected, the power spectrum of wing lift is much more 
sensitive to the lower gust frequency and decreases with increase in gust frequency. Note 
that the power spectrum is relatively insensitive to small b' values except at the larger 
frequencies . 

Figure 6 gives the percent reduction in the power spectra of figure 5 when the 
alleviation system is employed. A value of R = 100 indicates complete alleviation. 
Frequencies at which R = 0 on any of the curves indicate the point where the power 
associated with the alleviated wing first becomes larger than the power associated with 
the unalleviated wing. For any value of I*, there is a reduction in alleviation as gust 
frequency increases. This results because the high-frequency components of gust velocity 
over the wing span are not highly correlated with the velocity sensed by the vane. The 
correlation becomes worse as Z* increases. 

Figure 6 shows that the gust frequency range over which continuous alleviation is 
obtained (that is, R > 0) decreases with increase in wing size and vane location ahead of 
the wing. In all cases, however, the greatest reductions occur at the lower frequencies, 
which are more important for passenger comfort. 
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Lag Effects in Flap Response 

Notice that as the vane is moved farther ahead of the wing in figure 3 (that is, I* 
increases); the gust alleviation becomes less. The reason for this is that the flaps move 
too soon. If the vane is mounted ahead of the wing, there should be a time lag in the 
response of the flap to allow the gust sensed by the vane to reach the wing. This time 
lag can be incorporated into the system by varying the frequency of the servosystem 
driving the flaps. The present study always uses a damping ratio ? of 0.7, which means 
a response with small overshoot. 

Nondimensional me an- square lift .- The nondimensional mean-square lift $'g(0) is 
presented as a function of the nondimensional natural frequency k of the flap servo in 
figure 7 for b' values of 0.05, 0.1, and 1.0. This figure shows that there are optimum 
servo frequencies for minimizing the values of mean-square lift for given values of I* 
and b'. For example, if b' = 0.05, then the curve in figure 7(a) corresponding to 
I* = 0.5 shows that the optimum value of k for reducing the nondimensional mean- 
square lift is about 61. As k increases, the values of mean- square lift 5/g(0) approach 
the in-phase values shown in figure 3 for the particular values of I* and b' under 
consideration. 


•2r 


S 



(a) Id' = 0.05. 


Figure 


7.- Effect of servo response of the alleviation system on the 
mean-square value of e(x). £ =0.7. 
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Gain constant K .- Figure 8 gives the values of K used to generate figure 7. As 
k increases, the gain constant K approaches the in-phase values shown in figure 4 for 
the particular values of I* and b' under consideration. 

Optimum values of k and K .- Figure 9 presents the optimum values of servo fre- 
quency k and gain constant K required to minimize the nondimensional mean-square 
lift for different values of I* and b'. The resulting percentage reduction in the 

mean- square values of lift on the basic wing when these optimum values of k and K 
are used is shown in figure 10 . 

As an example in the use of figures 9 and 10, consider the following specific values: 
b = 15.24 meters (50 ft) 

U = 60.96 meters/sec (200 ft/sec) 

Z = 7.62 meters (25 ft) 

L = 304.8 meters (1000 ft) 



Nondimensional natural flap servo frequency, k 


(a) b' = 0.05. 

Figure 8 .- Gain constant used to minimize the mean-square value of e(x). 
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Figure 9-- Optimum values of K and 
k to minimize the mean-sq.uare 
value of e(x). 



Figure 10.- Percentage reduction in 
mean-square lift on the basic ving 
when optimum values of k and K 
are used in the flap servo system. 


Then, 


Z* = 0.5 


=0.05 


From figure 9, the optimum value of k for this value of I* and b' is about 61. From 
the definition of k, the dimensional natural frequency of the flap servosystem is 

related to k as 
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(On = 


k 



where the factor 2 tt is introduced so that cOn can be expressed in Hz (or cycles 
per second). 

Thus, the optimum natural frequency of the flap servo in the example is 


(0 


= 61 
" 27t(5) 


= 1.94 Hz 


The optimum gain constant K shown in figure 9 is very close to 1. Furthermore, 
the reduction in the mean-square lift, given by figure 10, is about 99 percent. 


A Specific Flap Control System 

By definition, the natural frequency coj^ corresponding to an optimum value of k 
depends on the unknown scale of turbulence L. In addition, the atmosphere is isotropic 
in patches at most and the scale of turbulence changes. Hence, it might be better not to 
try to optimize and K for a specific value of L, but rather to choose constant val- 
ues of ct)n K so that over a given range of possible L values, acceptable allevia- 

tion results will be obtained. Intuitively, it would appear that if such a system were fea- 
sible, the system would also function well in an atmosphere which is not quite isotropic. 

A fixed practical set of control parameters which were chosen somewhat arbitrarily 
for investigation are 

cijjj = 1.6 Hz 

K = 0.9 

? = 0.7 

Figure 11 shows the percentage reduction in the mean-square lift of the basic wing 
when this alleviation system is used for b' values of 0.05 and 0.1. The percentage 
reduction is presented as a function of the ratio L/U and vane location I*. The point 
of figure 11 is that the reduction is very substantial when the constant control parameters 
are used, and the unknown value of L/U falls within the wide range of values shown. 

For the airplane used in the flight tests of references 2 and 3, the wing span b was 
about 15.24 meters (50 feet) and the forward velocity U was about 60.96 meters/sec 
(200 ft/sec). Calculations were made to determine the reduction in the power spectrum 
of lift of the basic wing in this case. It was found that the results varied only slightly for 
reasonable values of L. The reason for this is that the wing span is much smaller than 
the scale of turbulence. The reduction in the power is shown in figure 12 as a function of 
gust frequency for different vane locations. The vane location for the flight tests (refs. 2 
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and 3) was about I = b/2 = 7.62 meters (25 ft). For this vane location, the power is 
reduced 100 percent at low frequencies and 70 percent at 2 Hz, which was considered to 
be the highest frequency of interest in the flight tests from the standpoint of passenger 
comfort. Since the alleviation system becomes less effective as the gust frequency 
increases, a low-band-pass filter should be used to attenuate the alleviation system 
beyond the frequency range of interest. Hence, it appears that if the frequency of the 
flap system is about 1.6 Hz, very effective gust alleviation can be realized for the type 
of aircraft used in the flight tests. 

In this study, the statistical properties of the lift on an airplane wing equipped with 
a gust alleviation system in flight through isotropic turbulence have been considered. A 
reduction in the wing lift due to the gusts by the alleviation system is, however, only part 
of the problem. A complete analysis should include the dynamic response of the entire 
airplane to the gust as was done, for example, in reference 1. Consideration of the span- 
wise variation of gusts has, however, indicated design parameters (flap servo gain and 
natural frequency and vane location) that give substantial reduction of the lift due to gusts. 

CONCLUDING REMARKS 

An analysis has been made of the effects of spanwise variations of gust velocity in 
isotropic turbulence on a gust alleviation system which employs an angle-of-attack vane 
mounted ahead of the wing to sense the vertical gust velocity. The wing flaps were moved 
in response to the vane deflection by a linear second-order servosystem to produce a lift 
opposite to that produced by the gust. Since the gust velocity actually varies across the 
wing span and the vane only senses the gust velocity at a single point, the flaps cannot be 
expected to eliminate the effects of the gusts completely by being deflected in response to 
a signal from the vane. The effectiveness of the alleviation system is examined by con- 
sidering the statistical properties of the residual lift on the alleviated wing. Considera- 
tion of spanwise variation of gusts has indicated design parameters (gain and natural fre- 
quency of the flap servosystem and vane location) that give substantial reduction of the 
lift due to gusts. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 11, 1971. 
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